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EXECUTIVE  SUMMARY 


CONTRACT  TITLE  AND  NUMBER; 

MHD  Compressor  -  Expander  Conversion  System  Integrated  with  a  OCR 
Inside  a  Deployable  Reflector. 

DE-AC22-88PC8B651 


CONTRACTOR  NAME  AND  ADDRESS: 

ANSALDO  S.p.A. 

Corso  F.M.  Perrone,  25 
16161  Genova  -  Italy 

START  DATE:  4-21-1988 

CDMMPLtTION  DATE;  4-20-1989 

This  work  originates  -from  the  proposal  "MHD  Compressor-Expander 
Conversion  System  Integrated  with  a  BCR  Inside  a  Deployable 
Re-flector"  submitted  by  ANSALDO  under  PRDA  RA-22-87PC90271 , 
January  1987.  The  proposal  concerned  an  innovative  concept  o-f 
nuclear,  closed-cycle  MHD  converter  for  power  generation  on 
space-based  systems  in  the  multi-megawatt  range.  The  basic 
element  of  this  converter  is  the  Power  Conversion  Unit  (PCU) 
consisting  of  a  gas  core  reactor  directly  coupled  to  an  MHD 
expansion  channel.  Integrated  with  the  PCU,  a  deployable 
reflector  provides  reactivity  control.  The  working  fluid  could 
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be  either  uranium  hexafluoride  or  a  mixture  of  uranium 
hexaf luoride  and  helium,  added  to  enhance  the  heat  transfer 

properties.  The  original  Statement  of  Work,  which  concerned  the 
whole  conversion  system,  was  subsequently  redirected  and 
focused  on  the  basic  mechanisms  of  neutronics,  reactivity 
control,  ionisation  and  electrical  conductivity  in  the  PCU. 
Furthermore,  the  study  was  required  to  be  “inherently  generic 
such  that  the  analysis  and  results  can  be  applied  to  various 
nuclear  reactor  and/or  MHD  channel  designs". 

Specific  goals  of  the  project  were: 

-  To  evaluate  the  performance  of  the  external  reflector. 

-  To  determine  the  influence  of  additional,  fission-induced 
mechanisms  on  the  ionization  levels  of  the  working  fluid. 

-  To  estimate  the  electrical  conductivity  levels  that  can  be 
attained  in  the  MHD  channel,  taking  into  account  the  effects  of 
F-  or  other  negative  ions. 

Two  major  conclusions  can  be  drawn  from  the  accomplishment  of 
this  project: 

-  A  reactor  configuration  has  been  obtained  which  avoids  to 
decouple  the  neutronics  of  the  cavity  and  its  internal 
reflector  from  the  external  shell,  thus  enabling  the  latter  to 
be  used  for  reactivity  control. 

In  particular  the  results  show  that  the  movable  reflector  is  able 
to  make  the  reactor  subcritical  by  1000  pern  at  least,  and  to 
control  the  reactor  from  CZP  to  HFP  condition,  accounting  for  a 
reactivity  change  up  to  3300  pem  associated  to  depletion. 
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-  The  conditions  expected  in  the  MHD  channel  show  that  the 
ionization  and  conductivity  levels  can  be  amplified  through 
'secondary'  processes  induced  by  the  'primary'  hot  electrons 
which  are  directly  generated  by  the  fission  reactions. 

The  numerical  calculations  suggest  that  conductivity  levels  in 
the  range  10-100  (fim)"1  are  consistent  with  the  character i sti cs 

of  short  nozzle,  high  power  density  generators.  In  fact,  the 
nozzle  length  turns  out  to  be  of  critical  importance  because  of 
the  sharp  decrease  of  residual  fission  and  gamma  densities  as  a 
function  of  the  distance  from  the  reactor  outlet. 

This  conclusion  remains  true  even  if  the  effect  of  dissociation 
and  attachment  are  included  in  the  numerical  model. 

Furthermore,  a  preliminary  evaluation  of  the  influence  of  wall 
material  and  thickness  shows  the  potential  for  the  fission  and 
gamma  densities  in  the  duct  to  be  increased  of  one  order  of 
maqnitude,  in  comparison  with  the  values  used  in  the 
conductivity  calculations,  if  a  material  with  favourable 
nuclear  properties  is  adopted. 
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1  INTRODUCTION 

i.i  Project  Objectives 

This  work  concerned  an  innovative  concept  of  nuclear, 
closed-cycle  MHD  converter  -for  power  generation  on  space-based 
systems  in  the  mui ti— megawatt  range.  A  schematic  diagram  of  the 
concept  is  shown  in  Fig.  1.1.  The  basic  element  of  this  converter 
is  the  Power  Conversion  Unit  (F'CU)  consisting  of  a  gas  core 
reactor  directly  coupled  to  an  MHD  expansion  channel.  The  working 
fluid  could  be  either  uranium  hexafluoride  or  a  mixture  of 
uranium  hexafluoride  and  helium,  added  to  enhance  the  heat 
transfer  properties. 

Two  major  items,  both  concerning  the  F'CU,  were  addressed  in  the 
project : 

1.  Because  of  the  high  levels  of  temperature  expected  in  the 
reactor  cavity,  reactivity  control  is  not  achieved  by  neutron 
absorbers  (control  rods)  located  in  the  cavity,  but  by 
moving  part  of  the  reflector,  located  outside  the  pressure 
vessel,  in  order  to  change  the  amount  of  neutrons  escaping 
from  the  core. 

2.  The  mere  levels  of  temperature,  however,  cannot  sustain 
levels  of  ionization  and,  hence,  of  conductivity,  capable  of 
providing  efficient  power  generation  in  the  MHD  channel. 
Fission  products  may  represent  a  possible  additional  source 
of  ionization. 
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This  report  documents  the  results  o-f  work  done  to  investigate 
these  aspects  of  the  orooosed  conceot.  ANSALDO  aodressed  both 
item  i.  ana  2.:  tne  latter  required  tne  scientific  support  of  the 
Universities  of  Florence  and  h'isa.  department  of  rhy^ics. 


1,2  Report  Organisation  . 

Section  3..  Reflector  Modeling,  illustrates  the  reactor 
conf i duration,  as  defined  throughout  contract  worn:,  Describes  the 
codes  and  procedures  adopted  in  the  analysis,  and  presents  the 
results  concernino  the  performance  of  the  external  movable 
reflector. 

Section  4. ,  working  Fluid  Characteristics,  investigates  tne 
f i ssi on- i noticed  ionisation  mechanisms  acting  in  the  working  gas, 
and  provides  a  system  of  simplified  rate  equations  for  numerical 
solution. 

Section  5.,  Electrical  Conductivity  Effects,  addresses  the  study 
of  the  electrical  conductivity  levels  in  the  MHD  channel ,  and 
presents  the  results  of  numerical  estimates  performed  with 
different  gas  and  duct  parameters. 

References,  Tables  and  Figures  are  located  at  the  end  of  each 
section . 
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2  PROJECT  DESCRIPTION 

This  section  is  based  on  the  contents  of  paragraph  4.2  in  the 
Statement  of  Work. 

The  work  is  divided  into  four  tasks. 

Task  1.  Reflector  Modeling 

This  task  shall  provide  the  nuclear  analysis  of  the  concept, 
performed  on  a  reference  PCU  geometry,  using  developed  and  open 
computer  codes  available  at  ANSALDG. 

The  reflector  modeling  shall  address  the  following  items: 

-  Reflector  reactivity  worth. 

-  Power  density  distribution  inside  the  PCU,  at  constant  core 
power  level. 

Task  2.  Working  Fluid  Characteristics 

This  task  shall  provide  an  analytical  model  enabling  the 
prediction  of  the  temperature  distribution  and  ionization  levels 
of  the  working  fluid  in  the  PCU. 

Since  the  equilibrium  thermal  ionization  is  not  expected  to 
maintain  sufficient  electrical  conductivity  for  efficient  power 
generation,  this  task  shall  investigate  additional  ionization 
mechanisms  that  can  be  present  in  a  fissioning  plasma. 
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Task  3.  Electrical  Conductivity  Effects 

This  task  shall  provide  appropriate  relations  enabling  the 
prediction  of  the  electrical  conductivity  of  the  working  gas.  For 
this  purpose,  this  task  shall  primarily  address  the  evaluation 
of  the  electron  mobility  in  the  gas.  The  possible  contribution 
(beneficial  and/or  detrimental)  to  electric  conductivity  of  the 
other  charged  particles  in  the  gas,  such  as  F-  ions,  shall  be 
taken  into  account. 

Task  4.  Reporting 

This  task  concerns  the  execution  of  the  reporting  requirements 
incorporated  in  the  contract. 
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3  REFLECTOR  MODELING 


3.i  Symbols,  Terms  and  Abbreviations 


a/ o 

BOL 

BU 

CZP 

f 

HFP 

HZF' 

<3 

pern 

Power  density 

? 

A? 


r 


Absorption 
Atom  per  cent 
Beginning  of  life 
Fuel  burn-up 
Cold  zero  power 
Fission 

Hot  -full  power 
Hot  zero  power 
Energy  group 

Percent  mille  (a  reactivity  change  o-f  1  pem 
equals  a  reactivity  change  of  10e-5  A^>) 

The  thermal  power  produced  per  unit  volume  of 
the  core  (W/cm*) 

Reactivity 

Change  in  reactivity,  defined  as 
=ln<k2/kl) ,  where  kl  and  k2  are  the 
eigenvalues  obtained  from  two  calculations 
that  differ  only  in  the  values  assigned  to 
the  independent  variables 
Removal 
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shutdown  margin  The  amount  o-f  negative  reactivity  (  ^>  )  by 

which  a  reactor  is  maintained  in 
subcri ti cal  state  at  CZF'  conditions  after  a 
control  trip 

Teff  Resonance  effective  temperature  of  the  fuel 

tr  Transport 
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3.2  Reference  Core  Description 

General 

The  results  presented  in  this  section  are  referred  to 
the  reactor  reference  conf iguration  shown  in  Figures  3. 1-3.2. 

This  conf iguration  is  a  simplified  arrangement  subject  to 
the  limitations  of  the  codes  adopted  in  the  analysis;  in 
particular,  a  coarse  discretization  of  both  composition  and 
geometry  has  been  made  at  the  core  inlet  and  outlet  conical 
segments. 

A  simple  axial-inlet  flow  pattern  has  been  adopted  in  the 
analysis;  however,  the  cooling  requirements  of  the  wall  made  it 
necessary  to  add  an  internal  diffuser  that  divides  the  inlet 
flow  into  two  parts,  the  larger  flowing  in  the  annulus  between 
the  internal  diffuser  and  the  external  one. 

In  this  reactor  the  same  fluid  acts  as  working  medium  and  fuel; 
it  flows  across  the  reactor  cavity  and  consists  of  highly 
enriched  gaseous  UF6  (region  1  of  Fig.  3.1):  90%  a/o  enriched 
in  U235  with  the  reactor  at  BOL. 

The  core  supplies  a  thermal  power  of  277  Mw  at  full  power. 

The  following  parameters  have  been  fixed  for  the  power 
density  distribution  calculations  inside  the  core  and  for 
the  reactivity  integral  worth  evaluation  of  the  movable 
reflector: 
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-  Reactor  at  full  power  (HFP) , 

-  BOL, 

-  no  Xenon,  no  5m, 

-  inlet  pressure  15  MPa, 

-  UF6  -flow  rate  1000  Kg/s, 

-  inlet  temperature  1720  K, 

The  rationale  -for  the  selection  of  these  system  parameters  is 
provided  in  section  5.2. 

Besides,  in  order  to  estimate  the  reactivity  requirements  during 
reactor  operations,  two  additional  core  conditions  have  been 
defined: 

a)  condition  I: 

-  Core  at  cold  zero  power  (CZP) , 

-  BOL, 

-  UF6  inlet  pressure  52  bar, 

-  UF6  inlet  temperature  573  K, 

-  UF6  flow  rate  1000  Kg/s, 

b)  condition  II: 

-  Core  at  hot  zero  power  (  HZP) , 

-  BOL, 

-  UF6  inlet  pressure  15  MPa, 

-  UF6  inlet  temperature  1720  K, 

-  UF6  flow  rate  1000  Kg/s, 

The  above  conditions  are  possible  if  the  reactor  control 


system  can  assure  that: 
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A)  condition  I  is  reached  without  using  nuclear  heating  to 
attain  the  gaseous  phase  of  the  -fuels 

B)  the  fuel  density  at  the  core  inlet  is  kept  constant  (0.36 
g/'cm~)  by  controlling  the  inlet  parameters  during  reactor 
operations. 

Components 

The  dimensions  of  the  reactor  components  (shown  in  Figures 
3.  1-3.2)  have  been  fixed  after  some  preliminary  considerations, 
based  on  the  thermal  and  mechanical  feasibility  of  the  system 
and  on  the  basic  neutronic  requirements  enabling  the  core 
reactivity  to  be  controlled  by  a  movable  reflector,  sliding 
axially  on  the  external  side  of  the  pressure  vessel. 

As  already  pointed  out  in  [13,  the  main  problem  is  the 
selection  of  a  pressure  vessel  material  highly  transparent  to 
neutrons. 

It  can  be  shown  that  the  adoption  of  the  most  common 
metallic  materials  ,  such  as  nickel  alloys  or  stainless-steel, 
introduces  a  strong  reactivity  penalty  and  decouples  the  inner 
shells  from  the  external  movable  reflector  to  such  a  degree 
that  reflector  worth  becomes  negligible. 

For  the  above  reasons  a  Zr/Nb  alloy  (UNS=R6090i>  has  been  finally 
selected  as  pressure  vessel  material. 


The  cylindrical  shell  thickness  (10  cm)  has  been  evaluated 

prel i minari 1 y  in  conformity  with  the  A5ME  III  Code  C21  for  a 

permanent  static  load  condition  at  15  MPa.  In  order  to 

minimize  creep  phenomena.  the  pressure  vessel  maximum 

temperature  has  been  supposed  to  be  lower  than  620  K. 

To  keep  the  pressure  vessel  temperature  at  such  a  low  value 
an  active  heat  removal  system  has  to  be  provided. 

Since  the  maximum  temperature  in  the  reactor  cavity  is  expected 
to  reach  about  2800  K  close  to  the  wall,  highly  refractory 
materials  have  been  used;  in  particular,  a  graphite  shell  (10  cm 
thick)  and  a  BeO  shell  (10  cm  thick)  have  been  interposed  between 
the  fuel  and  the  pressure  vessel. 

These  shells  have  to  be  cooled  by  helium  gas  flowing  across  a 
number  of  coolant  channels  obtained  in  the  shells  themselves. 

The  purpose  of  the  coolant  system  is  twofold: 

1)  to  remove  the  gamma  and  neutron  heating  from  the  inner  shells; 

2)  to  create  a  strong  temperature  gradient  between  fuel  and 
pressure  vessel.  Tab.  3.1  shows  the  radial  temperature  profile 
assumed  for  the  core  reference  configuration. 

Furthermore,  in  order  to  avoid  the  corrosion  reactions  between 
(JF6  and  graphite,  a  thin  Mo-alloy  liner  (0.6  mm  thick)  has  been 
placed  to  insulate  the  fuel  from  the  graphite  surface. 

An  adequate  liner  material  is  TZM  (UNS=F‘/M,R03640)  that  has  a 
high  melting  point  (3300  K)  and  a  liquidus  temperature  at 
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Typical  -fields  of  application  of  this  material  are  heat  engines, 
extrusion  dies  and  nuclear  reactors,  as  it  retains  good 
mechanical  performances  at  high  temperature. 

It  should  be  pointed  out  that  only  Mo-allovs  can  be  probably 
used  for  structural  parts  in  the  reactor  cavity  (e.g.  the 
diffuser);  however,  since  these  alloys  act  as  strong  absorbers 
in  the  proposed  concept  of  gas  reactor ,  their  use  should  be 
limited. 

In  this  preliminary  feasibility  analysis,  however,  the  internals 
(  diffuser  and  gas  coolant  tubes)  have  been  ignored;  the 
reactivity  penalty  due  to  the  internal  materials  can  be 
balanced  by  increasing  fuel  enrichment. 

Finally,  beryllium  has  been  selected  as  the  material  of  the 
movable  reflector. 

Mass  and  volume  of  the  reactor  components  are  summarized  in  Tab. 

*7  7 
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Analysis  Codes 

The  principal  computer  codes  adopted  in  the  analysis  are  GBC-4 
[33  (zero-dimensional ) ,  ANISN  [43  (one-dimensional),  SQUID-360 
[53  (two-dimensional),  WAPITI-GAS  [63  (two-dimensional), 
MERCURE-IV  [73  (three-dimensi onal ) .  A  brief  presentation  of  these 
codes  is  provided  below;  additional  information  can  be  found  in 
the  references. 
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6GC-4  Code: 

GGC-4  is  a  point  model  neutron  spectrum  program. 

The  fast  and  thermal  spectrum  calculations  are  performed  by 
the  methods  of  GAM  C 8 ]  and  GATHER  193,  respectively. 

Input  microscopic  cross  section  data  are  based  on  the 
ENDF/B-IV  Library  [103 

The  eneray  range  is  divided  into  a  number  of  fast  groups  (up 
to  99)  and  thermal  groups  (up  to  101). 

ANISN  Code  : 

ANISN  solves  the  one-dimensional  Boltzmann  transport  equation 
for  slabs,  and  for  cylindrical  or  spherical  geometries. 

As  secondary  calculation,  the  detailed  flux  generated  as 
the  solution  of  the  Boltzmann  equation  may  be  used  to  perform  a 
group  reduction  of  the  cross  sections. 

SQUID-360  Code: 

SQUID  is  a  multigroup  finite  difference  diffusion  code  in  two 
dimensions  (X-Y  or  R-Z). 

WAPITI-GAS  Code: 

WAPITI-GAS  is  a  finite  difference  code  developed  by  ANSALDO  for 
two-dimensional  (X-Y  or  R-Z)  thermal-fluid  analysis  of 
compressible  systems.  It  is  derived  from  the  French  code 
WAPITI. 
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MERCURE-IV  Code: 

The  Monte  Carlo  code  MERCURE-IV  computes  in  a  three-dimensional 
heterogeneous  geometry  heating  and  gamma  dose  rates,  and  by 
the  point-wise  kernel  attenuation  in  a  straight  line,  the  fast 
fluxes. 

Tabulated  accumulation  factors  give  the  contribution  of 
the  scattered  gamma  rays.  The  program  calculates  the 
accumulation  factors  for  mixtures;  for  multi-layer  media  the 
Kitazume  formulation  is  used. 

Additional  support  codes  have  been  used  for  special  calculations 
such  as  the  correlation  factors  and  the  cross  section  update 
during  the  iteration  process  between  SQUID  and  WAPITI. 


3.4  Features  of  Performed  Calculations 

Fuel  and  structural  material  cross  sections 

In  accordance  with  the  scheme  reported  in  Fig.  3.3  the  cross 
sections  have  been  calculated  by  using: 

a)  GGC-4  for  the  fuel  cross  sections  with: 

-  fast  energy  range  from  14.9  MeV  to  2.38  eV  and  B1 
approximation  (GAM); 

-  thermal  energy  range  from  2.38  eV  to  zero  and  BO 


O 


approximati on  (GATHER); 
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-  50  group  microscopic  cross  section  generation  (41  -fast 
groups  and  9  thermal  groups)  to  be  input  to  ANISN  code  in 
order  to  generate  -fuel  and  structural  cross  sections. 

The  assumed  energy  cuts  are  shown  in  Tab.  3.3. 

b)  GGC-4  as  in  item  a)  but  -for  the  structural  cross  sections 
with: 

-  B1  approximation  in  GATHER. 

c)  ANISN  code  using  cylindrical  geometry  and  F'0-54 
approximation  in  order  to: 

-  investigate  the  -fuel  5  group  cross  section  dependence  on  the 
temperature  and  density  of  the  fuel; 

-  generate  the  structural  material  macroscopic  cross  sections 
and  a  5  group  cross  section  data  set  to  calculate  the  fuel 
macroscopic  cross  sections  to  be  input  to  SQUID  code  as 
a  function  of  temperature  and  density  at  each  iteration  step. 

The  assumed  energy  cuts  in  the  range  14.9  MeV  -  0.0  eV  are 
14.9  MeV,  0.821  MeV,  150  KeV,  5.53  KeV,  0.625  eV,  0.0  eV. 

Reactor  Calculations 

In  accordance  with  the  scheme  reported  in  Fig.  3.4,  total 
reflector  worth  and  the  distributions  of  thermal  power, 
fuel  temperature,  and  fuel  density  have  been  calculated  by  using: 
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a)  SQUID-360  code  with: 

-  simplified  R-Z  geometry  and  material  compositions, 

-  k-eff  calculation, 

-  spatial  -fluxes  calculations  in  order  to  evaluate  the  thermal 
power  density  distribution  to  be  input  to  WAPITI  code. 

b)  WAPITI-GAS  code  -for  the  thermal-fluid  analysis  with: 

-  R-Z  geometry, 

-  turbulent  flow, 

-  density  and  temperature  distributions  of  the  fuel  to  be  input 
to  SQUID-360  Code. 


.;.,5  Nuclear  Cross  Section  Calculations 


The 

macroscopic  cross  sections  required  as 

input 

data  for 

the 

full 

core  R-Z  SQUID  calculations  have 

been 

generated 

by 

using 

the  procedure  outlined  in  Fig.  3.3. 

The 

cross  sections 

are 

referred  to  : 

-fuel 

? 

O 


-structural  material. 
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Fuel  cross  sections 

It  has  been  assumed  that  inside  the  reference  core  configuration 
(see  par.  3.1)  fuel  temperature  and  density  may  change  in  the 
range  1800-3000  K  and  0.15-0.45  g/cm3,  respectively ,  as  shown  in 
Fig.  3.5. 

A  sensitivity  analysis  has  been  performed  by  applying  the 
procedure  of  Fig.  3.3  to  each  of  the  points  1-9  of  Fig.  3.5. 

The  results  are  summarized  below: 

a)  the  fuel  microscopic  cross  sections  change  only: 

-  with  fuel  density  in  the  4th  energetic  group  (epithermal 
range) , 

-  with  fuel  temperature  in  the  5th  energetic  group  (thermal 
group) ; 

b)  the  5  group  structural  material  cross  sections  are  not 

significantly  affected  by  the  fuel  conditions  under  the 
assumption  that  the  radial  temperature  profile  in  the 

structural  material  doesn't  change. 

Consequently,  in  order  to  update  the  fuel  macroscopic 
cross  sections  during  the  iterative  SQUID  calculations  of 
reactor  power  di stribution  with  thermal-fluid  feedback,  the 
following  correlations  can  be  used: 


g  g  n 

E  (T,n)  =  Z  (To, no)  — 
j  j  no 


g=l ,3  and  j=tr,a,f 
or 

g=l,4  and  j=r 
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4  4  n  4 

Z  (T,n)  =  Z  (To, no) —  +AZ  (n)  j=tr,a,f 

j  j  no  j 

5  5  n  5 

Z  (T,n)  =  Z  (To.no)—-  +AZ  (T)  j=tr,a,f 

j  j  no  j 


4 

AZ  (n)  =  ntajn^  +  bjn  +  cj)  j=tr,a,f 

j 


5 

AZ  (T)  =  n(AjTr  +  BjT  +  Cj) 
j 


where: 


[n]=g/cm-'  E  i  ]=K 


no=0. 3 


g/cm-2' 


To=2l00  K 


9 


The  correlation  coefficients  aj ,  bj,  cj,  A j ,  B j ,  Cj  have  been 
estimated  by  using  the  previous  ANISN  calculations. 

The  results  are  reported  in  Tab.  3.4. 
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Structural  material  cross  sections 

The  structural  macroscopic  cross  sections  Tor  the  full  core 
R-Z  SQUID  calculations  have  been  evaluated  by  ANISN  with  the 
procedure  shown  in  Fig.  3-3  by  fixing  fuel  density  at  O.oo  g/cm- 
and  fuel  temperature  at  2100  K,  which  represent  average  values 
in  the  core  reference  conf iguration. 


3.6  Reactor  Calculations 

The  reactor  calculations  have  been  performed  by  applying  the 
scheme  reported  in  Fig.  3.4  to  the  core  reference  conf iguration 
with: 

a)  external  reflector  out,  65%  of  the  fuel  flowing  in  the 
annular  diffuser  and  35%  of  the  fuel  flowing  in  the  internal 
diffuser; 

b)  as  in  item  a),  but  with  the  external  reflector  100%  in; 

c)  external  reflector  100%  in,  70%  of  the  fuel  flowing  in  the 
annular  diffuser  and  30%  of  the  fuel  flowing  in  the  internal 
diffuser. 


The  proposed 

scheme  can  be 

useful 

to 

study  a 

generic 

stationary  core 

conf iguration  in 

order 

to 

evaluate 

power 

density  distribution  and  k-effective,  along  with  fuel  temperature 
and  density  distributions. 

In  short,  the  calculation  procedure  is  the  following: 

1)  a  guess  power  density  distribution  is  fixed, 
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2)  WAF'ITI-GAS  computes  the  -fuel  density  and  temperature 
distributions  generated  by  the  input  power  density 
distribution, 

3)  the  fuel  density  and  temperature  distributions  are  used  to 
update  the  nuclear  macroscopic  cross  sections  for  each 
fuel  composition, 

4)  SQUID  computes  a  new  power  density  distribution  and 
eigenvalue, 

5)  a  test  on  the  convergence  of  the  main  variables  is  made: 
if  it  is  successf ul ,  than  the  iterative  process  is  stopped; 
if  not,  it  returns  to  step  2). 

Figures  3. 6-3. 8  report  the  axial  profile  of  power  density  at 
different  radii  for  the  above  configurations. 

Tables  3. 5-3. 7  show  the  respective  temperature  distributions  in 
the  reactor  cavity.  The  values  are  computed  at  the  center  of  each 
mesh.  The  axial  coordinate  is  given  in  the  first  column  of  the 
table,  the  value  at  the  top  corresponding  to  the  reactor  inlet: 
radially,  the  meshes  are  numbered  from  the  reactor  axis. 


3.7  Reactivity  Requirements 

In  this  core  concept,  excess  reactivity  is  controlled  by  the 
external  movable  reflector.  This  excess  reactivity  includes: 
a)  reactivity  changes  due  to  changes  in  gas  working  temperature 
between  the  cold  and  hot,  no  load  conditions: 
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b)  reactivity  changes  due  to  changes  in  reactor  power  over  the 
power  range  of  operations; 

c)  reactivity  changes  due  to  fuel  depletion; 

d)  reactivity  associated  to  the  minimum  shutdown  margin. 

The  reactivity  change  associated  to  each  item  is  reported  in  Tab. 
3.B. 

The  integral  worth  of  the  movable  reflector  has  been  evaluated 
by  means  of  the  reactor  calculations;  an  uncertainty  of  about  107. 
has  been  associated  to  the  results. 

The  reactivity  change  between  CZF'  and  HZP  does  not  appear 
to  be  significant  because  in  this  reactor  the  Doppler  effect  has 
been  considered  only  for  U236,  as  in  the  GAM  library  there 
isn't  any  information  about  the  Doppler  effect  in  U235. 

Finally,  since  it  is  not  in  the  scope  of  this  preliminary 
analysis  to  make  accurate  assumptions  about  the  system  operating 
modes,  it  has  been  supposed  that  the  reactor  is  operated 
at  full  power  level  continuosly. 

It  has  been  estimated  that,  with  a  3300  pern  of  reactivity  change 
available  for  depletion,  reactor  life  is  limited  to  a  few 
hours  because  of  the  strong  Xenon  poisoning  that  affects  this 
type  of  core.  This  time  is  a  function  of  the  UF6  overall 
plant/core  mass  ratio  as  shown  in  Fig.  3.9. 

The  full  power  level  can  be  restored  after  an  interval  of  time 
whose  length  depends  on  the  operational  strategy;  the  process  can 
be  repeated  until  fuel  depletion  allows  it. 
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Fig  3.10  shows  the  evolution  of  reactivity  after  a  shutdown 
performed  at  the  end  of  the  first  lifetime  at  full  power,  whose 
length  is  about  5  hours.  In  this  case,  which  is  referred  to  a 
UF6  mass  ratio  of  5,  the  time  required  for  reactivity  to  become 
positive  again  is  approximately  30  hours  from  shutdown:  this 
waiting  time  is  reduced  by  increasing  the  mass  ratio. 

The  minimum  value  of  reactivity  that  occurs  9  hours  after 
shutdown  in  Fig.  3.10,  corresponding  to  a  maximum  of  the  Xe 
concentration ,  is  explained  by  the  presence  of  a  source  of  Xel35, 
the  radioactive  decay  of  1135,  in  competition  with  its 
disappearance  due  to  decay  to  Csl35. 
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Tab.  3.1  Radial  temperature  profile  in  the  core 
reference  configuration. 


Radial  Range  (*) 

(cm  ) 

Material 

Average  Tempe 

(K) 

50.0-53.5 

Graphite 

2200 

53.5-56.5 

II 

1500 

56.5-60.0 

II 

1200<T<1500 

60.0-63.5 

BeO 

900<T<1200 

63.5-66.0 

II 

900 

66.0-70.0 

•1 

600<T<900 

70.0-80.0 

Zr-Nb 

<600 

80.0-110.0 

Be 

<600 
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Tab. 

3 . 2  Mass 

and  volume  summary 

for  the 

core 

reference  configuration. 

Reactor  Material 

Average 

Volume 

Mass 

%  of 

Density 

(dm3) 

(kg) 

Total  Mass 

(kg/dm3 ) 

UF6 

0.29 

1723 

500 

1.9 

T2M 

10. 1G 

2.93 

20 

O.l 

Graphite 

1.71 

1320 

2257 

0.5 

UcO 

3.02 

1317 

3977 

14 . 0 

Zr-Nb  alloy 

6.5 

1 G  3  4 

10G21 

30.8 

Dc 

1.85 

5012 

9272 

34.8 

Total  Reactor  Mass 


2G65G  100.0 
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Tab.  3.8  Reactivity  balance  for  the  core 
reference  configuration. 

Movable  reflector  integral  worth 

Assumed  uncertainty  (10%) 

From  HZP  to  HFP 

From  CZP  to  HZP 

Shutdown  margin  plus  depletion 


9.6% 


-0.96% 


-4.3% 


0% 

4.34% 
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n  ( g / cm3 ) 


T :  temperature 
n:  density 


Fig.  3.5  Fuel  characterization  as  a  function  of  temperature 
and  density. 
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Fig.  3.7  Power  density  vs.  axial  < 
coordinates  . 
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Fig.  3.10  Reactivity  vs.  time  (UF6  mass  ratio  =  5). 
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4  WORKING  FLUID  CHARACTERISTICS 


4.1  Symbols,  Terms  and  Abbreviations 


n 

G 

v 

No 

a 

T 

ne 

ni 

no 

rei 

8 


k 

m 

A 

M 

e 

r 


eo 

Tp 


Number  of  particles  per  unit  volume 
Collision  cross  section 
Particle  velocity 

Number  of  neutral  molecules  per  unit  volume 

Recombination  coefficient 

Temperature 

Number  of  electrons  per  unit  volume 

Number  of  ions  per  unit  volume 

Number  of  primary  electrons  per  unit  volume 

Thermal i ration  time  between  electrons  and 

i  on  s 

Quantity  in  the  expression  for  thermal i ration 
time 

Boltzmann’s  constant 
Electron  mass 
Ion  mass  in  amu 
Ion  mass 
Electron  charge 

Quantity  in  the  expression  for  coulomb 
collisions 

Energy  of  primary  electrons 
Temperature  of  primary  electrons 
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Ts 

Ti 

ree 

Tpei 

rsei 

I 

R 

Np 

a,b 

IS 

Vg 


Temperature  of  secondary  electrons 
Temperature  of  ions 

Thermal ization  time  between  primary  and 
secondary  el ectrons 

Thermal ization  time  between  primary  electrons 
and  ions 

Thermal i z at i on  time  between  seco.noary 

electrons  and  ions 
Ionization  potential 
Ionization  degree 

Number  of  gas  molecules  per  unit  volume 

Coefficients  in  the  rate  equations 

Ratio  of  no  to  Nq 

Gas  ve 1 oc i t  y 

Axial  coordinate 

Average 
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4.2  Phase  I 

This  phase  provided  a  qualitative  assessment  of  the  types  of 
primary  reactions  considered  as  the  most  important,  and  the 
evaluation  of  the  relevant  cross-sections,  using  current 
literature  data  and  theoretical  models. 

The  Production  of  Primary  Electrons 

Whenever  235U  undergoes  a  fission  reaction  a  number  of  ionizing 
agents  are  generated.  Fission  fragments  are  the  most  important  as 
they  absorb  around  80*/.  of  the  energy  released.  The  energy  balance 
of  a  fission  reaction  is  illustrated  in  TABLE  1  (data  are 
essentially  similar  for  thermal  or  fast  fission  reactions  and  for 
the  two  uranium  isotopes) Cl]. 

TABLE  1 


1.  Fission  Fragments  165  MeV 

2.  Prompt  gamma  rays  (bremmstrahlung)  7  MeV 

3.  Delayed  gamma  rays  (discrete  spectrum)  6  MeV 

4.  <i+  and  (3-  rays  7  MeV 

5.  Neutrons  5  MeV 

6.  Neutrinos  10  MeV 

TOTAL  200  MeV 
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Excepting  neutrinos,  that  are  not  stopped  in  the  gas  and  escape 
■from  the  core  because  of  their  extremely  small  cross-section, 
all  the  other  remaining  products  are  stopped  within  a  few  mms  of 
gas  (the  neutrons  travel  the  furthest  distance  but  are,  however, 
partly  recovered  by  the  reflector,  and  the  gamma  rays  occur  as 
electron-positron  pairs  in  a  few  cms)  and,  therefore,  they  lose 
their  energy  partly  by  ionizing  the  constituent  elements  of  the 
system  and  partly  through  the  fission  reactions. 

In  this  phase  we  examined  the  elastic  and  inelastic 
cross-sections  of  the  various  fission  products  on  gaseous  UF6 
and  He  C23.  Objective  of  the  analysis  was  the  evaluation  of  the 
number  of  ionizations  that  take  place  for  each  fission, 

neqiecting  the  effects  induced  by  this  so  called  primary 
ionization. 

The  enerpy  produced  by  the  fission  processes  is  largely  converted 
into  el ectrons: the  majority  of  these  (equivalent  to  about  165 
MeV)  are  due  to  ionization  processes  induced  by  the  fission 
fragments,  and  have  a  diffused  spectrum  with  an  average  energy  of 
500  eV. 

A  smaller  part  (20  MeV)  results  from  gamma  ray  conversion  and 
beta  decays  and  may  have  an  energy  of  iMeV  per  electron.  However , 
in  this  case,  the  energy  is  rapidly  distributed  over  other 
electrons,  and  hence  it  can  be  assumed,  to  a  good  approximation, 
that  a  single  fission  actually  produces  4. 10e5  primary  electrons 
characterized  by  an  average  energy  of  500  eV. 
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These  primary  ionizations  are  mostly  made  up  by  "stripping"  an 
electron  without  modifying  the  molecular  or  atomic  structure  of 
the  target. 


4.5  Phase  II 

This  phase  addressed  the  study  of  the  rate  equations  in  the 
reactor  core.  providing  an  approximate  analytical  solution,  and 
identifying  the  critical  parameters. 

The  Production  of  Secondary  Electrons 

The  primary  electrons  give  rise  to  processes  such  as  elastic  and 
inelastic  scattering  (excitation  of  bonding  electrons  in  the 
molecule),  secondary  ionization,  and  recombination. 

Elastic  scattering  is  low  at  energies  greater  than  10 
eV  and  can  be  almost  neglected  in  most  cases.  The 
inelastic  scattering,  on  the  other  hand,  is  a  very  important 


reaction 

channel  since 

it  absorbs  around  2/3 

of 

the  primary 

electron 

energy.  This 

value  is  difficult  to 

be 

computed  or 

measured.  It  depends  on  the  energy  and  on  the  structure  of  the 
molecules;  data  available  in  literature  give  a  ratio  between 
ionization  and  inelastic  collision  cross  sections  of  0.2  -  0.5 
[33.  These  reactions  are  part  of  the  gas  thermal ization  and 

heating  process. 
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The  atomic  cross-section  for  ionization  processes  on  UF6  around 
500  eV  is  of  the  order  of  10e-15  cm* .  The  cross-sections  are 
known  with  errors  of  25'A.  It  should  be  noted  that  the  ionization 
cross-sections  on  He  have  a  very  high  threshold  value,  due  to  the 
high  ionization  potential  (about  24  eV)  ,  and  a  very  low  value, 
due  to  the  closed  shell  structure  of  He  which  permits  elastic 
scatterina  to  predominate.  Hence,  the  addition  of  He  is  expected 
to  increase  the  electrical  conductivity  to  a  very  small  extent. 
The  contribution  to  the  ionization  process  is  given  by  the 
product  of  the  number  of  electrons,  the  cross-section,  and  the 
velocity.  This  yields,  therefore: 

dn/dt  =  +nGKv>Ng 

where  Q  is  the  cross  section,  numerically  equal  to 
10e-15,  <v>  is  the  velocity  in  cm/sec,  equal  to  <v>=6. 10e7TT 

(where  T  is  the  temperature  given  in  eV)  and  Ng  is  the  number  of 
molecules  per  cm3. 

The  probability  of  ionization  by  collision  should  be  actually 
obtained  by  averaging  the  cross-section ,  which  depends  upon  the 
electronic  energy,  over  the  assumed  Maxwell  velocity  distribution 
[43.  However,  the  numerical  value  indicated  for  Q  is 
correct  for  the  order-of-magni tude  estimate  of  this  phase. 
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Recombination 

The  neutral  ization  probability  of  the  ionized  system  is 
related  to  the  product  of  the  electronic  density,  the 
ionic  density,  and  the  density  of  a  third  element  essential 
for  energy  and  momentum  conservation  [43.  This  third  element  may 
be  either  a  neutral  atom  or  an  ion,  or  even  an  electron.  The 
higher  mobility  of  the  electrons  makes  the  following  the 
preferred  recombination  process  in  our  case  : 

ne  +  ne  +  ni  — >  ne  +  ni 

This  parameter  is  extremely  sensitive  to  the  density  and 
temperature  of  the  system. 

The  rate  of  recombi nati on  can  be  expressed  by  the  following 
relation: 

dn/dt  =  -  ocn2 

where  a  =  10e-25n/TA- E,  n  is  the  electron  (or  ion)  density  per 
cm*  and  T  is  the  temperature  expressed  in  eV. 

Note  the  dependence  upon  n*  (the  coefficient  oc  is  linearly 
related  to  n) ,  which  reminds  that  recombi nation  is  a  three-body 
process  [51. 
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The  dependence  on  T  is  very  strong  because  the  electron  velocity 
is  the  other  key  parameter  to  the  problem,  both  when  calculating 
the  cross-section  and  when  successively  calculating  the 
average  on  the  velocity  distribution,  necessary  to  determine  the 
probability  per  second  of  the  process. 

Simplified  Rate  Equation 

The  effects  of  recombinati on  and  ionization  may  be  combined  into 
a  rate  equation  for  the  electronic  part  of  the  system: 

dn/dt  =  nQ<v>Ng  -  an2 

Substituting  the  numerical  values  into  the  above  equation 
yields: 


dn/dt  =  n  3. 10el3  TJT  -  10e-25  n?7T4' c 

Ng  was  assumed  to  be  5. 10e20  molecules  per  cm*. 

Let  us  now  suppose  to  be  under  steady  state  conditions 
(hence  neglecting  dn/dt). 

This  assumption  requires  that  a  unique  equilibrium  temperature 
exists,  at  least  among  the  various  ionized  species,  and  that  the 
time  to  reach  such  conditions  is  much  shorter  than  the  time 
during  which  the  gas  remains  in  the  reactor  core.  In  fact, 
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electrons  around  500  eV,  electrons  of  a  few  eVs  (secondary) ,  and 
ions  at  almost  the  same  temperature  of  the  gas  coexist  during  a 
first  phase. 

However,  as  we  shall  see  better  afterwards,  the  conditions  in  the 
reactor  core  lead  to  extremely  short  thermal ization  times 
compared  with  the  transit  time  of  the  gas,  and  hence,  in  spite  of 
the  simple  approach,  this  first  evaluation  may  be  considered 
quite  realistic. 

Let  us  now  introduce  the  energy  balance,  that  is  determined  by 
the  distribution  of  the  ionization  and  thermalization  energies: 

ynoSOO  =  n(T+I)  =>  n(T+10)  =  y2.S-10e21 

where  y<l  is  the  part  of  energy  involved  in  ionization  processes 
and  I  is  assumed  to  be  equal  to  10  eV,  a  value  typical  of  most 
atoms  and  molecules. 

The  rate  equation,  once  inserted  into  the  balance  equation, 
gives: 


To-'=.  (T+10)  =  165y  =>  T  =  2.  By-® 

Similarly,  the  implicit  value  of  the  electronic  density 
can  be  found: 


n  =  6-10e5(2.8-10e20y-n)°^ 
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which  yields  approximately 


n  = 

2.B- 10e20y 

It  should  be 

noted 

that 

the  recombination  process  is 

active, 

above 

all , 

at 

low  temperatures.  However,  the 

process 

temperature 

is 

determined  critically 

through 

the 

parameter 

y. 

The 

conventional  thermal 

region 

(T<4000  C)  is  indicated  by  values  of  y  less  than  10e-2. 

Therefore,  the  critical  parameter  that  must  be 
calculated  is  the  actual  distribution  of  energy  between  the 
ionization  and  thermal ization  processes  or,  if  one  prefers,  the 
electron  equilibrium  temperature. 

Thermal izati on  Times 

What  remains  to  be  seen  is  whether  the  secondary  electrons,  which 
still  have  sufficiently  high  energy,  of  the  order  of  2-3  eV,  can 
thermal ize  (and  hence  recombine)  rapidly  or  not. 

The  electron-ion  thermalization  time  decreases  with  the 
temperature  according  to  the  following  relationship  [53: 

e 

rei  =  — —  (kT)3',= 

n 


where  is  given  by  [53  as: 
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i 


3  m 1 /=  1632  A 

6  =  i"  . . 

8 (2tt)  e*  logT 

where  m  is  the  electron  mass,  e  the  electron  charge,  A  the  ion 
mass  in  a.m.u.,  and  logT  a  parameter  which  varies  slightly  with 
density  and  temperature  (it  is  related  to  the  smallest  collision 
parameter  in  a  Coulomb  interaction) .  In  the  case  of  UF6,  the 
value  of  9  is  3.10e9,  with  kT  expressed  in  eV.  Assuming  that  the 
•final  energy  has  a  value  of  a  few  eVs,  with  a  density  of  about 
10el9  el/cm3,  the  electron-ion  thermalization  time  is  of  the 
order  of  10e-9  s. 

Another  important  parameter  is  the  electron-electron 
thermalization  rate.  In  fact,  in  our  case,  high  energy  electrons 
(around  500  eV)  and  secondary  electrons  (around  10  eV)  will 
coexist  for  a  certain  period  of  time.  The  electron-electron 
thermal i zati on  time  is  reduced  by  the  electron/ion  mass  ratio, 
that  is  ree=(l/1832A)rei .  Both  thermal ization  times  are  much 
shorter  than  the  gas  replacement  time  (0.1  s,  typically)  and, 
hence,  it  is  reasonable  to  assume  steady  state  conditions  to 
calculate  the  concentrations  of  the  chemical  species  in  the 
reactor. 
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Complete  Rate  Equations 

The  "complete"  rate  equations  take  actually  into  account  only  a 
tew  of  the  numerous  possible  reactions  in  the  system.  On  the 
basis  ot  the  cross-section  experimental  data  and  ot  the  analysis 
reported  in  the  previous  paragraphs,  we  selected  the  reactions 
which  were  considered  as  the  most  significant  for  the 
purpose  of  this  study. 

reactions  such  as  the  multiple  ionization  of  atoms  or  molecules 
have  been  neglected.  In  other  words,  only  UFi,  UFi+  (with  i=0,6> 
and  F-  type  systems  are  considered  and  possible  states  such  as 
UFi-,  UFi++,  etc.  are  neglected. 

Under  this  assumption,  there  are  16  unknown  parameters  in  the 
problem: 

ni  =UFi/cmT-  (  7  variables  ), 
ni+=UFi+/'cmT-  (  7  variables  ), 
nf  =  F/cm3, 
nf2  =  F2/cmT', 

nf-  =  F-/cm3  (3  variables), 

in  addition  to  n  that  represents,  as  before,  the  number  of  free 
electrons.  Therefore  16  equations  are  required  to  solve  the 
system  as  a  function  of  three  basic  input  parameters:  Ng  (gas 
density  in  the  reactor) ,  no  (number  of  primary  electrons)  and  the 
average  primary  radiation  energy  eo  =  500  eV. 

The  rate  equations  provide  17  relations,  one  for  each  variable: 
the  energy  balance  yields  the  missing  condition  (system 
neutrality  is  implicit  in  the  rate  equations). 
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Method  of  Solution 

Be-fore  proceeding,  it  is  necessary  to  discuss  some  assumptions 
that  should  be  introduced  in  a  numerical  approach  to  the 
evaluation  o f  ionization  levels  and  gas  composition  in  the 
cavity. 

1)  The  system  is  homogeneous. 

This  is  a  good  approximation  given  the  parameters  of 
the  reactor  and  the  properties  of  the  system  being  examined. 
In  fact,  the  mass  transfer  rates  are  very  small  in  comparison 
with  two  significant  fluid  parameters,  the  thermal 
transfer  rate  and  the  speed  of  sound?  in  addition,  the 
thermal i ration  processes  take  place  on  small  space 
scales,  of  the  order  of  a  few  cms.  Hence,  there  is  no  use  in 
introducing  the  convection  terms  into  the  rate  equations,  and 
considering  the  unknowns  as  a  function  of  position. 

2)  The  system  is  under  steady  state  conditions. 

This  assumption  is  undoubtedly  true  when  considering  the 
thermal ization  time  scale  (10e-9  s) ,  and  allows  the  rate 
equations  to  be  simplified  considerably  by  neglecting  the  time 
dependence  of  the  variables. 

3)  Cross-section  values. 

In  our  approach,  there  are  46  cross-sections  (which  are  strongly 
energy  dependent)  constituting  a  set  of  input  parameters  to  the 
problem.  They  are  not  variables  but  physical  quantities  that  can 
be  experimentally  measured.  In  any  case,  these  parameters  are  not 
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well  known  and  it  is  therefore  necessary  to  create  a  simple 
theoretical  model  to  predict  the  cross-section  relative  weights 
and  their  dependence  on  energy. 

The  system  can  be  solved  by  means  of  successive  iterations.  A 
tentative  value  is  assumed  for  the  number  of  electrons  (for 
example  the  value  calculated  with  the  simplified  equation  shown 
previously)  and  the  system  solved  with  this  value.  The  solution 
in  this  case  is  very  simple  because  the  system  is  almost  linear 
since  the  majority  of  the  quadratic  functions  are  related  to  the 
presence  of  the  unknown  n.  Now  the  average  energy  is  calculated, 
the  values  of  the  cross-sections  are  redefined  and  the  number  of 
electrons  recalculated  using  the  first  rate  equation.  The  value 
thus  obtained  is  inserted  back  into  the  equations  and  the 
procedure  repeated.  After  a  few  iterations  the  result  begins  to 
converge  towards  the  exact  value. 


4.4  Phase  III 

This  phase  addressed  the  study  of  the  rate  equations  in  the  MHD 
channel,  providing  a  simplified  set  of  equations  for  numerical 
solution. 

The  following  assumptions  have  been  introduced  to  describe  the 
evolution  of  ionization  degree  and  temperature  in  the  MHD 
channel : 


00123 


1)  gas  density  is  uni-form  during  the  transit  time  in  the 
channel ; 

2)  a  uniform  distribution  of  fission  reactions  provides  a 
uniform  background  of  primary  electrons  (500  eV) ; 

3)  as  a  first  approach,  the  secondary  species  are  supposed  to 
derive  essentially  from  the  ionization  of  UF6.  Therefore,  the 
maximum  level  of  charge  density  cannot  overcome  gas  density. 
The  consequences  of  this  statement  will  be  discussed  iater; 

4)  the  loss  of  electrons  is  ascribed  to  three-body 
recombi nati on : 

5)  since  the  local  neutrality  condition  is  well  verified 
because  of  the  very  short  Debye  length  compared  with  plasma 
dimensions,  it  is  assumed  ne=ni  for  secondary  charges. 

Even  a  simplified  approach  requires  to  solve  the  system  giving 
the  thermal  exchange  rates  among  the  ionized  species.  Since 
charge  density  is  variable,  the  energy  balance  must  be  coupled  to 
the  rate  equation  for  electron  production: 


dTp 

no  . — 

dt 


noTp  -  neTs  noTp  -  neTi 


ree 


rpei 


(1) 
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ali 
ne  — 

at 


naip  -  neu  neis  -  neTi 

__________  +  i  .I... . . 

rpei  rsei 


(3) 


arte  a  exp(-I/Tp)  exp(-I/Ts) 

_  =  -  No(  - no  +  - - ne)  -  bne-Ts”*' c  (4) 

dt  I  Tp1'^  Ts1''^ 

Equations  (1),  (2)  and  (3)  represent  the  cooling  ot  primary 
electrons  by  tnermai  exchange  with  secondary  electrons  and  ions, 
the  neatino  ot  secondary  electrons  by  exchange  with  the  primary 
ones  and  the  transfer  ot  energy  -from  secondary  electrons  to 
ions  (Eq. (2) ) ,  and  finally  the  heating  of  ions.  Id,  Ts  ana 
Ti  are  the  temperatures  (in  eV>  of  primary  electrons,  secondary 
electrons,  and  ions,  respectively.  no  is  the  uniform  ana 
constant  density  of  5bu  eV  electrons,  ne=ni  is  the  charge 
density  produced  by  additional  ionization. 

rpei,  rsei, and  ree  represent  the  thermalization  time  constants 
for  the  involved  interactions.  The  general  form  of  r  has  been 
diven  in  the  previous  paragraph.  Therefore,  the  time  constants  of 
Equations  (1),  (2)  and  (3)  are  given  by: 


e 

ree  =  —  (kTp)--'r  and  rpei  =  1632ree 
ne 
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6 

rsei  =  1832A  -  (kTs)3'= 

ne 

In  Eq.  (4)  ,  the  -first  term  represents  the  rate  of 
ionization  produced  by  primary  and  secondary  electrons, 
and  the  second  one  the  three-body  recombination  o-f 
secondary  charges,  the  -fissions  providing  the  gas  with  a 
continuos  source  of  primary  electrons. 

I  is  the  ionization  potential,  a  and  b  are  numerical 
coefficients,  and  No  is  the  density  of  neutral  molecules,  related 
to  ne  by  the  mass  conservation  relation. 

It  is  useful  to  introduce  the  degree  of  ionization  as  a 
variable,  namely  R-ne/Ng  where  Ng  is  the  gas  density.  The  density 
of  neutral  molecules  is  replaced  by  No=Ng-ne,  and  the  density 
of  primary  electrons  by  no=flNgj  no  is  given  by  the  number  of 
f i ssi ons/cm3s  times  the  ionization  efficiency  of  each  fission, 
whose  value  has  been  estimated  in  a  previous  paragraph  at  about 
4.10e5.  Therefore,  the  system  of  Equations  < 1 )  —  < 4 )  is  transformed 
in 


dTp  RNg  R 

dt  9Tp-=  0 

dTs  0Ng  R 

dt  ©Tp'-'2  ft 


m  Tp  R 

-  _( - >*'=  (Ts-Ti)  -] 

M  Ts  0 

m  Tp  R 

-< - '>*'=  (Ts-Ti)  -I 

M  Ts  0 


(5) 


(6) 
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dTi  m  r^Ng  R  m  Tp  R 

_  =  -  -  [Tp  -  -  Ts  -  -( - )*'  =  (Ts-Ti)-]  (7) 

dt  M  eTp~'=  (5  M  Ts 


dR  aNg(l-R)  exp(-I/Tp>  exp<-I/Ts)  bNg=K~ 

__  =  [  j %  +  . . . R]  -  .  (8) 


at 


I 


I  5  : 


IS 


The  assumption  of  a  uniform  density  distribution,  related 
to  the  rather  flat  profile  of  gas  density  during  the  transit 
in  the  channel,  allows  one  to  neglect  the  contribution  of  the 
terms  Vg.grad(T)  and  Vg.grad(n)  in  comparison  with  the  time 
derivatives.  Therefore,  an  elementary  volume  of  fluid  is 
modified  in  its  properties  before  appreciable  deplacement 
takes  place.  As  a  consequence,  the  integration  in  time  is 
equivalent  to  follow  the  space  evolution,  with  the  variable 
change  x=Vg.t.  It  turns  out  that  a  non-uniform  source  of 
primary  electrons  can  be  included  in  the  model,  under  the 
condition  that  the  induced  gradients  still  verify  the  inequality 
Vg.  grad  (T)  «dT/dt  ana  Vg.grad  (n)<<dn/dt.  For  instance,  in  the 
presence  of  a  step  distribution  of  the  fission  source  in 
a  length  xo,  it  is  sufficient  to  introduce,  in  the  solution 
of  the  system,  the  condition  no=0  for  t>xo/Vg.  Alternati vely , 
the  system  can  be  solved  in  the  space  variable  by 
introducing  d:<=Vg.dt.  In  general,  it  is  not  strictly 
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required  that  gas  velocity  and  density  are  uniform  and 
constant:  it  is  sufficient  that  significant  variations  of  Vg 
and  density  occur  in  a  time  longer  than  the  time  required  to 
reach  a  stationary  state  for  Tp,  Ts,  Ti,  and  R. 
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5  ELECTRICAL  CONDUCTIVITY  EFFECTS 


5,1  Symbols,  Terms  and  Abbreviations 


tr 

Me 

F'g 

?C 

J 

L 

D 

u 

B 

Ms 

? 

st 

T 1 
0 

0o 

d 

x 

Ue 

Ui 

E 


Electrical  conductivity 

Mass  o f  the  superconducting  coil 

Power  generated  by  the  MHD  channel 

Conductor  mean  density 

Permeability  of  free  space 

Current  density 

Channel  length 

Channel  diameter 

Gas  velocity 

Magnetic  field  strength 

Mass  of  the  magnet  structure 

Density  of  the  structural  material 

Working  stress  level  of  the  structural 

material 

Gas  temperature  at  the  compressor  inlet 

Neutron  flux  along  the  duct 

Neutron  flux  at  the  entrance  of  the  duct 

Diameter  of  the  entrance  section  of  the  duct 

Axial  coordinate  along  the  duct 

Electron  velocity 

Ion  velocity 

Total  electric  field 
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Eo 

Eint 

e 

m 

M 

Vet 

Vit 

k 

Te 

Ti 

ne 

ni 

n 

Q 

Ng 

R 

e 

Di 

De 

pe 

pi 

Ts 

c 

QN 

Tp 


External  electric  field 
Internal  electric  field 
Electron  charge 
Electron  mass 
Ion  mass 

Total  collision  frequency  for  electrons 
Total  collision  frequency  for  ions 
Boltzmann's  constant 
Electron  temperature 
Ion  temperature 

Number  of  electrons  per  unit  volume 
Number  of  ions  per  unit  volume 
Number  of  particles  per  unit  volume 
Collision  cross  section 
Number  of  gas  molecules  per  unit  volume 
Ionization  degree 

Quantity  in  the  expression  for  thermal ization 
time 

Diffusion  coefficient  of  electrons 
Diffusion  coefficient  of  ions 
Electron  mobility 
Ion  mobility 

Temperature  of  secondary  electrons 
Coefficient  in  the  conductivity  relation 
Collision  cross  section  with  neutrals 
Temperature  of  primary  electrons 
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a,b 


no 


§ 

nf 

t 

T 

I- 


Temperatu.re  of  ions 
Coefficients  in  the  rate  equations 
Number  of  primary  electrons  per  unit  volume 
Ratio  of  no  to  Ng 

Number  of  primary  electrons  per  fission 
Number  of  fissions  per  unit  volume 
Time 

Time  constant 
Attachment  energy 


a- 


Coefficient  in  the  attachment  equation 
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5.2  Channel  Characteristics  and  Conductivity 

To  proceed  further  in  the  analysis,  it  is  necessary  to  combine 
conductivity  calculations  with  actual  channel  (and  nozzle) 
characteristics.  In  fact,  the  results  of  the  former  depend 
directly  on  the  latter,  as  channel  and  nozzle  lengths,  along  with 
the  density  distribution  of  the  gas  flowing  through  them,  are 
required  to  solve  the  rate  equations  (l)-<4)  or  (5)-(8)  given  in 
section  4.4.  Furthermore,  the  same  characteristics  influence  the 
distributions  of  residual  fission  neutrons  and  gamma  rays 
available  in  the  duct  which,  in  turn,  provide  the  source  of 
primary  electrons  to  be  introduced  in  the  rate  equations.  On  the 
other  hand,  the  channel  operating  characteristics  are  largely 
dependent  on  the  conductivity  levels  assumed  in  the  analysis.  For 
this  reason,  a  parametric  approach  has  been  adopted,  and  the 
generator  performance  has  been  estimated  for  two  levels  of 
conductivity,  namely  <r=10  and  tr=loO  dim)-1.  The  nozzle  and 
channel  conf igurations  arising  from  this  process  have  been  tested 
to  verify  if  they  were  consistent  with  the  conductivity  levels 
resulting  from  the  rate  equations.  Channel  performance  has  been 
evaluated  under  the  following  assumptions: 

-  Faraday  configuration  of  the  generator  electrodes  with  a  load 
factor  of  0.9; 

-  Straight  channel  with  a  constant  aperture  angle  of  7°; 

-  Nozzle  inlet  temperature  and  pressure  of  2300  K  and  15  MPa, 


respectively; 
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-  Channel  -flow  rate  of  1000  Kg/s; 

A  1-D  thermal -f 1 ui d  model  of  the  channel  and  nozzle  ducts  has 
been  aoooted  in  the  analysis.  For  eacn  of  the  two  reference 
levels  of  conductivity,  a  number  of  cnannei  conf igurations  have 
been  compared  as  a  first  approach  by  using  the  specific  mass 
relation  given  by  Rosa  in  [1]  for  a  channel  of  length  L  and 
diameter  D: 

Me  16  ^>c  L  1  1 

Pg  jjo  J  D  Pg  tto-lFB1 

where  yc  is  the  conductor  mean  density  and  J  is  the  current 
density.  This  relation  gives  the  mass  of  a  superconducting  coil 
per  unit  of  Generated  power.  In  the  conditions  of  the  analysis, 
tne  mass  of  the  structure  reouired  to  hold  the  coil  together 
should  be  comparable  with  the  coil  mass  given  by  Eq.  (1),  A  lower 
limit  for  the  former  can  be  estimated  with  the  following 
relation: 


f  f 

Ms  =  -  j— ^dv  (2) 

stj  JJO 


where  p  is  the  density  and  st  is  the  working  stress  level  of  the 
structural  material.  This  relation  too  is  derived  from  Cl  3;  the 
structure  unit  mass  for  a  saddle  coil  is  expected  to  be  within  a 
factor  2  of  the  value  given  by  Eq.  (2). 

In  addition  to  the  two  levels  of  conductivity  mentioned  above, 
the  conf igurations  that  have  been  investigated  are  featured  by 
different  values  of  the  inlet  Mach  number  and  induction  field  of 
the  generator  duct.  The  parameters  assumed  in  the  analysis  enable 
to  extract  power  in  the  120-150  MWe  range,  limiting  the  induction 
field  strength  at  reasonable  values  (6  T  at  most),  observing  the 
constraints  on  the  length  to  diameter  ratio  imposed  by  heat  and 
pressure  losses,  and  obtaining  a  specific  mass  (from  Eq.  (1))  in 
the  range  of  10e-2fl0e-3  Kg/KWe.  This  range  is  in  the  reach  of 
both  reference  levels  of  conduct! vity;  however,  a  decrease  of  tr 
must  be  offset  by  a  corresponding  variation  of  other  terms  in  Eq. 
(1),  In  particular,  the  values  of  u  and  B  are  assumed  as 
independent  parameters,  while  L/D  and  F'g  are  derived  from  the 
calculations.  It  is  evident  that  a  higher  level  of  conductivity 
makes  it  possible  to  attain  the  desired  range  of  power  density 
with  lower  values  of  B  and,  especially,  of  u  (and,  hence,  of  the 
Mach  number  at  the  generator  inlet).  This  latter  condition 
appears  to  be  of  special  importance,  as  a  lower  gas  velocity  in 
the  channel  means  a  shorter  nozzle  before  the  channel.  On  the 
other  hand,  the  distance  between  the  generator  and  the  reactor 
cavity  has  a  strong  influence  on  the  actual  level  of 
conductivity.  The  results  show  that,  if  one  pursues  a  high  power 
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density  condition  (that  is,  essentially,  a  high  value  of  <ru2> , 
the  system  should  work  at  a  high  conductivity,  low  velocity 
condition,  rather  than  the  opposite. 

It  should  be  pointed  out  that  such  influence  is  not  connected 
with  the  effect  of  the  field  on  the  reactor  core,  which  is 
negligible,  due  to  the  different  order  of  magnitude  of  the 
fission  ion  gyro  radius  in  comparison  with  their  mean  free  path. 
In  fact,  the  nozzle  length  turns  out  to  be  of  critical  importance 
because  of  the  sharp  decrease  of  residual  fission  and  gamma 
densities  as  a  function  of  the  distance  from  the  reactor  outlet. 
Finally,  a  few  considerations  about  the  parameters  adopted  in  the 
analysis. 

The  value  of  nozzle  inlet  temperature  (2300  K>  has  been  partly 
assumed  on  the  basis  of  material  requirements,  by  analogy  with 
conventional  MHD  converters,  and  partly  to  maintain  an  adequate 
level  of  cycle  efficiency  with  a  high  heat  rejection  temperature 
(T 1=600-900  K,  T1  being  the  compressor  inlet  temperature) ,  as  a 
space-based  system  is  expected  to  require.  Within  this  range  of 
T1  and  with  the  other  cycle  parameters  indicated  in  this  section, 
efficiency  may  vary  between  0.12  and  0.36.  It  should  be  stressed 
that  cycle  efficiency  is  heavily  influenced  by  heat  rejection 
temperature  and  regenerator  performance;  both  characteristics 
should  be  ultimately  selected  on  the  basis  of  minimum  system 
mass,  rather  than  maximum  efficiency. 
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The  nozzle  inlet  pressure  (15  MPa)  is  equal  to  the  average 
operating  pressure  of  the  reactor,  which  was  derived  on  the 
basis  of  an  average  reactor  temperature  and  density  of  2000  K  and 
0.3  g/cm3,  respectively.  The  latter  is  the  result  of  a 
preliminary  estimate  of  criticity.  The  same  calculations  gave  the 
cavity  diameter  (1  m)  assumed  in  the  neutron  analysis. 

The  flow  rate  of  1000  kg/s  results  from  the  power  output  to  be 
generated  by  the  converter,  assumed  in  the  range  120-150  MWe 
qtoss  and  20—30  MWe  net.  The  latter  is  actually  a  lower  limit 
applicable  to  the  highest  rejection  temperature  assumed  in  the 
calculations;  the  net  output  increases  to  60-80  MW  when  Tl  is 
lowered  to  600  K. 

The  temperature  of  the  gas  at  the  reactor  inlet  (1770  K)  has  been 
obtained  through  a  preliminary  analysis  of  the  working  cycle,  and 
it  is  typical  of  regenerative  Brayton  cycles  operating  between 
the  temperature  levels  mentioned  above.  The  same  analysis  gave 
the  reference  value  of  277  MW  for  the  thermal  power  provided  by 
the  reactor. 


5.3  Fission  and  Gamma  Density  in  the  Channel 

The  fission  and  gamma  power  densities  in  the  MHD  expansion 
channel  must  be  estimated  to  solve  the  rate  equations  and 
evaluate  the  conductivity  levels. 


MD  00123 


Fission  Power  Density 

The  -fission  power  density  calculations  have  been  carried  out  oy 
means  of  the  following  approximations: 

1.  Preliminary  evaluation  with  a  simplified  method  (see  Ref,  [21, 
pag,  290). 

2.  SQUID-  360  C31  diffusion  calculations  with  a  simplified 

two-dimensional  R-Z  geometry  and  walls  fabricated  from  A1203 
or  BeO. 

In  order  to  estimate  the  power  density,  the  duct  has  been 

divided  into  a  finite  number  of  regions  with  different  UF6 

density  (see  Fig.  5.1  and  Fig.  5.2). 

For  each  region,  the  fission  macroscopic  cross  sections  have 
been  derived  from  the  core  cross  sections  calculated  during 
Task  1.  reduced  by  the  appropriate  density  factors. 

The  neutron  flux  level  inside  the  MHD  channel  has  been 
calculated  as  transmission  through  a  cylindrical  duct. 

Preliminary  evaluation.  The  neutron  fluxes  at  the  inlet  of  the 
nozzle  have  been  derived  from  the  nuclear  analysis  of  Task  1. 

The  effect  of  wail  material  reflections  and  working  fluid 

absorptions  have  been  neglected. 

The  uncoil ided  fluxes  along  the  duct  have  been  calculated 
with  the  formula: 

d* 

0  =  0o  — — 

8x2 
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where  0o  is  the  -flux  at  the  entrance  of  the  nozzle,  d  is  the 
diameter  o-f  the  entrance  section  o-f  the  nozzle,  and  L  is  the 
distance  between  the  entrance  o-f  the  nozzle  and  the  point  where 
the  flux  is  calculated. 

The  results  <MeV/cm3s>  are  reported  in  Tab.  5.1. 

Diffusion  Calculations.  In  order  to  estimate  the  influence  of 
wall  material  and  thickness  on  fission  power  density,  a  second 
analysis  of  the  MHD  channel  has  been  carried  out  by  the 
SQUID-360  diffusion  code,  using  the  R-Z  geometry  to  model  the 
Power  Conversion  Unit  (gas  reactor  +  nozzle  and  MHD  channel),  and 
the  cross  sections  derived  from  Task  1.. 

The  conf  i  gu.rati  on  extends  radially  to  include  the  movable 
reflector  and  axially  from  the  core  midplane  to  1200  cm  in  the 
channel  (see  Fig.  5.2). 

As  to  the  spatial  discretization,  114  axial  and  37  radial 
meshes  have  been  used. 

The  duct  has  been  divided  into  9  regions  with  different  UF6 
densities  carried  out  from  Fig.  5.3. 

The  nozzle  and  channel  walls  have  been  assumed  to  be  fabricated 
from  A12D3  or  BeO.  For  each  material,  two  values  of  wall 
thickness  have  been  examined:  3  cm  and  6  cm.  The  adoption  of  a 
two-dimensional  model  doesn't  account  for  the  presence  of  the 
electrodes.  However,  it  shouldn't  be  unreasonable  to  assume  an 
electrode  wall  consisting  of  "thin"  electrodes  mounted  on  "thick" 
insulating  blocks;  furthermore,  the  electrodes  themselves  might 
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be  -fabricated  -from  zirconia  based  materials,  and  zirconium 
possesses  favourable  properties  being  a  low  capture  cross-section 
material,  which  would  allow  to  exploit  the  potential  of  a 
moderating  layer. 

As  previously  mentioned,  the  calculations  have  been  carried 
out  by  using  the  cross  section  data  from  Task  1.  collapsed 
to  a  5  group  structure. 

The  fission  power  for  each  region  is  shown  in  Fig.  5.4.  The 
results  of  the  cases  examined  are  reported  in  Tab.  5.2. 

Gamma  Power  Density 

The  gamma  fluxes  in  the  MHD  channel  have  been  calculated  by 
the  MERCURE  IV  point  kernel  code  [41. 

The  gamma  source,  consisting  of  fission  gamma,  is  located  in  the 
core. 

Build-up  factors  have  been  used  in  the  MERCURE  calculation. 

The  geometry  of  the  analyzed  conf iguration  is  shown  in 

Fig. 5.1. 

The  core  has  been  divided  into  3  regions  and  the  duct  into  6 
regions  with  different  UF6  density. 

Region— wise  UF6  densities  in  the  duct  have  been  carried  out 
from  Fig.  5.3. 

The  nozzle  and  channel  walls  are  assumed  to  be  fabricated  from 
A1203  material,  but  they  have  not  been  taken  into  account  for 
streaming  calculations. 
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The  calculations  have  been  performed  by  using  the  BIP  G2 
library  cross  section  data  from  Ref.  C51. 

The  results  are  reported  in  Table  5.3. 


5.4  Conductivity  Calculations 

A  Relation  for  Conductivity 

Owing  to  the  range  of  gas  transit  time  in  the  channel  <~  10 
ms)  and  to  the  very  small  Debye  length  compared  with  the 
dimensions  of  plasma,  the  motion  of  both  charged  species  can  be 
described  by  neglecting  the  convective  derivatives  of  the 
velocities  and  assuming  the  neutrality  condition  ne=ni. 
Therefore,  the  velocities  Ue  and  Ui  are  given  by  C6] 

eE  kTe  grad(n) 

Ue  - - - - <1> 

mVet  mVet  n 


eE 

MVit 


kTi  grad(n) 


MVit  n 


(2) 


£2100  aw 
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where  E  represents  the  total  electric  -field,  given  by 
the  superposition  of  an  external  field  Eo  and  the  internal  field 
Eint  qenerated  by  the  relative  motion  of  charges,  Vjt  is  the 
total  collision  frequency  including  momentum-exchange 
collisions  with  neutrals  and  between  electrons  and  ions,  namely 

RNg 

Vjt  =  <QUj>Ng(l~R)  +  -  (3) 

0Tes'= 

where  the  symbols  have  been  previously  defined.  The  internal 
field  is  obtained  by  equating  the  fluxes  of  ions  and 
electrons,  in  accordance  with  the  neutrality  condition.  It  turns 
out  that  in  the  absence  of  an  external  field,  the  current  is 
zero.  The  internal  field  is  given  by 

Di-Be  grad(n) 

Eint  =  — . . — 

pe+pi  n 

where  Dj=kTj/mjVjt  is  the  diffusion  coefficient  and  pj=e/mjVjt 
is  the  mobility.  Therefore,  current  density  becomes 


e=n 

J  =  en(Ui-Ue)=  e~n (pe+pi )Eo  =  -  Eo  (4) 

mVet 
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where  the  mobility  of  ions  has  been  neglected,  because  of  the 
very  small  mass  ratio.  By  using  the  same  notations,  the 
conductivity  is  given  by 

e“RTs -■/:r 


mCcTs--  ( 1  — R )  +R/h] 


where  Ts  is  in  eV  and  c  is  proportional  to  the  cross  section  with 
neutrals.  In  MKS  units 


lOe-ZRTs’^ 


0-  z  . ■ —  — —  (fim)"1  (6) 

3. 10e7SNTs=(l-R)+10e-4R 

QN  is  the  geometrical  cross  section  with  neutrals  in  ctrF.  If  the 
contribution  of  collisions  with  neutrals  is  negligible,  the 
conductivity  is  independent  on  the  degree  of  ionization  R  and 
increases  as  Ts  increases.  With  a  typical  value  of  QN-10e-15, 
the  conductivity  is  only  function  of  Ts  if  the  condition 


3. 10e-4Ts= 

R  »  - - 

1+3. 10e-4Ts= 
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is  verified.  In  case  of  total  ionization  <R=1),  the 
conductivity  is  simply  equal  to  10e2Ts3/=.  If  Ts£l0  eV,  a  degree 
of  ionization  of  107.  is  sufficient  to  give  <r=3. 10e3  (ftm)-1. 

Numerical  Results  and  Discussion 

The  system  of  Equations  (5) - (S)  derived  in  section  4.4  has  been 
solved  by  the  Runge-Kutta  fourth  order  method.  The  initial 
conditions  are  Tp=500  eV,  Ts=l  eV,  Ti=0.5  eV  and  R=o.  They 
correspond  to  typical  mean  values;  however,  the  stationary 
values  don't  depend  on  the  initial  conditions.  The  a  and  b 
parameters  in  Eq. (8)  are  chosen  on  the  basis  of  data  available 
in  literature,  namely  a=10e-6  and  b=10e-22  (cgs  units).  The 
value  of  b  is  probably  more  realistic  than  that  given  previously: 
in  any  case,  an  overestimation  of  the  recombination  rate  gives 
the  lowest  limit  of  achievable  ionization.  The  ionization 
potenzial  of  UF6  is  £15  eV.  The  thermal ization  coefficient 
is  10e4  (cgs  units).  The  initial  rate  of  primary  electrons  ft 
is  related  to  the  number  of  fissions  i.e  i3=§nf/Ng, 
where  §  represents  the  number  of  500  eV  electrons  per  fission, 
whose  value  has  been  estimated  at  4. 10e5.  Figures  5. 5-5. 7 
show  the  time  evolution  of  temperatures  and  degree  of 
ionization  at  different  initial  gas  densities.  It  has  been 
assumed  i3=10e-2,  that  corresponds  to  nf  =2. 5el0-8Ng. 
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Under  the  conditions  of  Fig. 5. 5  (Ng-5. I0el7  molecules/cm-) , 
the  system  reaches  a  stationary  conf iguration  in  less  than  1 
ps,  with  Tp=24  eV,  Ts=Ti=3.7  eV  and  R=0.067.  The  resulting 
conductivity  is  (r=7.10e2  (ftm)-1. 

Figures  5.6  and  5.7  show  the  system  evolution  with  Ng=3. lOelS 
molecules/cms  and  Ng=10e20  molecules/cm3. 

The  stationary  data  referring  to  Fig.  5.6.  are  Tp=140  eV,  Ts=Ti=9 
eV  and  R=.15,  with  nf=7.5.10el0  fissions,  which  leads  to 
(T=3.10e3  (i'im)"1,  whereas  in  the  case  of  Fig. 5. 7.  Tp=50  eV, 

Ts=T i =7  eV  and  R=0.02,  but  at  t=50  ns  (nf =2. 5. 10el2) . 

The  results  suggest  that  the  increase  of  neutral  density 
cannot  be  always  considered  as  a  factor  that  inhibits  the 

growth  of  ionisation  rate.  Fig.  5.8  shows  the  stationary  values 
of  the  relevant  parameters  as  a  function  of  the  initial  gas 
density. 

Fig.  5.8  shows  that  the  highest  values  of  conductivity  are 

obtained  in  a  region  around  10el8  molecules/cm3.  Of  course,  this 
value  depends  on  the  choice  of  the  input  parameters. 
Nevertheless,  a  variation  of  one  order  of  magnitude  of  the 
ionization  and  recombination  coefficients  (a  and  b)  slightly 
affects  the  results. 

The  production  rate  of  primary  electrons,  due  to  the 
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fission  reactions,  plays  a  relevant  role.  If  an  initial  rate  ft 
of  10e-3  and  an  initial  density  of  3. 10el8  molecules/cm3  are 
assumed,  the  stationary  values  become  Tp=21  eV,  Ts=Ti=2.4  eV  and 
R=0. 01 ,  i.e.  tr  =2. 10e2  (ftm)-1  instead  of  3. 10e3,  if  the  initial 
value  of  ft  is  10e-2. 

The  distribution  of  fission  rate  in  the  channel  is  very 
important  in  order  to  maintain  suitable  conditions  for  MHD 
conversion.  In  fact,  if  the  ionization  source  stops  at  a.  time 
to  (equivalent  to  a  sharp  cut  in  the  spatial  distribution  at 
xo=Vg. to),  the  further  evolution  of  the  ionization  rate  R  is 
given  by 


dR  bNa-R-- 


dt  Ts''-'3 

the  secondary  electrons  being  too  cold  to  sustain  the  ionization 
level.  This  equation  can  be  integrated  by  assuming  as  initial 
values  those  determined,  in  the  stationary  state,  by  the  fission 
source.  Thus,  the  time  evolution  is  represented  by 
R=Ro/  ( 1+t/r)  1 y'=,  where  r=Ts'T''=/2b  (RoNg)  =. 

By  introducing  the  data  of  Fig.  5.6,  R  is  reduced  by  a  factor  lu 
in  0.05  ps.  It  results  that  in  the  absence  of  a  ionization 
source,  the  conductivity  drops  quickly  to  unsuitable  values. 

The  time  behavior  of  temperatures  and  ionization  rate 
greatly  simplifies  the  analysis  of  the  fuel  chemical  composition, 
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with  special  regard  to  the 

formation  of 

negative  ions  and 

to  the 

contribution  to  ionization 

of  the 

species  arising 

from  LIF6 

dissociation.  Dissociation  is  important  from  the  point  of  view 
of  the  effective  neutral  density.  Reminding  the  results  shown 
in  Fig.  5.8,  an  increase  of  neutral  density  due  to  dissociation 
does  not  affect  dramatically  the  relevant  parameters.  On  the 
contrary,  if  the  initial  density  is  5. 10el7  molecules/cm3'-,  as 
it  could  be  expected  during  the  gas  expansion  in  the  MHD  channel, 
a  growth  of  a  factor  7,  due  to  dissociation  and  subsequent 
ionisation,  could  lead  to  an  improvement  of  the  final  degree 
of  ionization  and  temperatures,  whose  values  should 
approximate  those  of  Ng=3.5.10ei8  molecules/cm3  (see  Fig.  5.8). 
However,  accordina  to  the  measurements  of  Compton  C73,  the 
thresholds  for  the  formation  of  UFn+  type  ions  by  electron 
impact  ionization  are  respectively  14  eV  for  UF5+,  18  eV  for 
UF4+,  22  eV  for  UF3+,  26  eV  for  UF2+,  32  eV  for  UF+  and  40  eV  for 
U+.  Since  the  stationary  temperature  of  secondary  electrons  and 
ions  is  in  the  range  of  5  eV,  relevant  contributions  can  be 
only  expected  from  UF6  and  UF5. 

The  influence  of  the  formation  of  negative  ions  on  the 
conductivity  can  be  estimated  by  simple  arguments.  The  affinity 
potential  is  around  5  eV  with  a  cross  section  of  3  C71 
(the  behavior  is  similar  for  UF6-  and  F-).  By  averaging  over  a 
Maxwell  distribution,  the  contribution  of  attachment  in 
Eg. (8)  can  be  expressed  as 
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a- 

<QV>ne (NF+NUF) =ne (NF+NUF)  -  exp <-I-/kTs)  (cm’s)-1 

Ts1-'- 

where  a-  slOe-7  (cgs  units)  and  I— as— 5  eV.  By  introducing  this 
term  in  Eq. (8) ,  it  turns  out  that  tne  previous  results 
are  sliahtly  modified.  even  in  case  of  total  dissociation, 
i.e.  NF+NUF=6Ng.  In  particular,  the  degree  of  ionization 
decreases  of  less  than  107.  in  comparison  with  the  value 
obtained  without  attachment. 

From  the  above  discussion,  it  seems  that  the  model  of 
Equations  (5)-<8)  gives  a  reliable  description  of  the 
ionization  in  conditions  of  relatively  low  gas  density.  The 
model  does  not  consider  inelastic  collisions  between  electrons 
and  neutral  atoms,  with  excitation  of  internal  levels,  and 
transfer  of  enerov  between  ionized  and  neutral  gas.  The 
influence  of  these  effects  can  be  estimated  by  comparing  the 
collision  frequency  with  neutrals  and  between  electrons  and 
ions.  With  the  previous  notation,  the  loss  of  energy  by 
collisions  with  neutrals  can  be  neglected  if  (see  Eq. (6)) 

R 

-  »  3. lOel lQNTs= 

(1-R) 
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With  QN=10e-15  cmr,  Ts=9  eV,  R»2.5. 10e-2.  The  value  of  R=0. 15 
(Fig.  5.6)  ensures  that  only  about  10  7.  of  energy  is  tranferred 
to  the  neutral  gas.  As  a  consequence,  the  chemical 
composition  of  the  neutral  gas  will  be  simply  determined  by  the 
thermal  equilibrium  relations,  calculated  at  the  temperature  of 
the  neutral  gas,  being  negligible  the  heating  due  to  ionization. 
Fia.  5.9  represents  the  chemical  composition  of  UF6  calculated 
in  conditions  of  thermal  equilibrium  at  a  pressure  of  10  atm 
(from  Ref . C6 3 ) . 

The  percentaae  of  neutrals  of  differents  species  can  be 
directly  obtained  from  Fig.  5.9  at  the  temperature  of  the  gas 
flowing  in  the  channel,  whose  value  is  not  modified  by 
the  ionization  processes.  Referring  for  instance  to  the 
conditions  of  Fig.  5.6,  the  stationary  state  in  the  channel  can 
be  approximately  estimated  as: 

907.  neutral  gas,  of  composition:  95  7.  UF6,  5  7.  UF5  (at  Tg=20OU 
K) . 

107.  ionized  yields,  of  composition:  507.  UF6+,  507  UF5+,  F-<10  7. 
(from  the  experimental  results  of  Ref.C73) 

With  an  initial  density  of  3.  lOelB  molecules/cm"',  we  obtain: 

2.565. 10el8  molecu.les/cm3  of  UF6,  0.  135.  lOelB  molecuies/cm3  of 
UF5,  2.7. 10el7  molecules/cm3,  1.5.i0el7  molecules/cm3  of  UF6+  and 
UF5+,  less  than  3. 10el6  molecules/cm3  of  F-. 

Finally,  Figures  5.10-5.12  report  the  results  obtained  with  the 
gas  density  distribution  of  Fig.  5.3,  and  the  fission  density 
distribution  of  Tab.  5.1.  The  latter  is  conservati ve:  the 
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preliminary  evaluation  of  the  influence  of  wall  material  and 
thickness  reported  in  Tab.  5.2  shows  the  potential  for  the 
fission  and  gamma  densities  in  the  duct  to  be  increased  of  one 
order  of  magnitude  at  least,  in  comparison  with  the  values  used 
for  the  conductivity  calculations,  if  a  material  with  favourable 
nuclear  properties  is  adopted. 
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'nab.  5.1  Fission  power  density  in  the  MHD 
channel,  simplified  calculations. 


Distance  from 
Entrance  of 
Nozzle  (  cm  ) 


Fission  Power  Density 
(  MeV  /  cm3  sec  ) 


50 

1 

1 

I 

1.72 

E+l  2 

200 

1 

1 

1 

2.90 

E+10 

400 

1 

1 

| 

2.80 

E+09 

600 

1 

1 

1 

9.20 

E+08 
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Tab.  5.2  Fission  power  density  in  the  MHD 
channel,  diffusion  calculations 
(MeV/cm3  s) . 


A1203 
Thickness 
3  cm  I 


BeO 

Thickness 
c  j  3  cm  |  d 


b  6  cm 


s 

Aver . Power 
Dens i ty 

Ratio 

a/a 

Ave  r . Power 
Density 

Rat io | 
b/a  | 

Aver . Power 
Dens  i  ty 

|  Rat io 

1  c/a 

Aver . Power 
Density 

|  Ratio 

1  d/a 

1 

1 . 10E+14 

1.0 

1 . 35E+14 

1.22  | 

2.10E+14 

|  1.9 

3.20E+14 

1 

|  2.9 
| 

2 

6 . 54E+12 

1.0 

9.61E+12 

1  .47  | 

2 . 10E+13 

i  3.2 

4.70E+13 

1 

|  7.2 

I 

3 

9 . 45E+11 

1.0 

1 . 85E+12 

1 . 96  | 

3.00E+12 

|  3.2 

1.20E+13 

1 

| 12 . 76 

4 

1 .93E+11 

1.0 

4.62E+11 

2.4  j 

7 . 51e+ll 

|  3.9 

4.08E+12 

1 

|  21.2 

I 

5 

4.42E+10 

1.0 

1.23E+11 

2.781 

2 .08E+11 

j  4.7 

1 . 40E+12 

1 

|  31.7 

I 

6 

7.99E+09 

1.0 

2. 56E+10 

3.2  | 

4.95E+10 

j  6.2 

4.08E+11 

j  51.1 

I 

7 

4 . 25E+08 

1.0 

1.70E+09 

4.0  j 

3.91E+09 

1  9.2 

4.19E+10 

1 

|  98.6 

1 

8 

1.83E+08 

1.0 

8.22E+08 

4.5  j 

1.96E+09 

|  10.7 

2.36E+10 

1 

|  129.2 

I 

9 

1 .37E+08 

1.0 

6.55E+08 

4.8  | 

1 .58E+09 

111. 6 

1 .96E+10 

1 

|  143.3 
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Tab.  5.3  Gamma  power  density  in  the 
MHD  channel . 


I  Distance  from  |  Gamma  Power  Density 

|  Entrance  of  I 

|  Nozzle  (  cm  )  |  {  MeV  /  cm3  sec  ) 

i  _ i _ 
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Graphite 
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BeO 
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Zr/Nb  Alloy 
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Be 
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7 

8 
9 


14 


Fig.  5.1  Three  dimensional  model  for 
fission  density  simplified 
calculations  and  Mercure 
calculations. 
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TZM 

2 

12 

Graphite 

3 

13 

BeO 

4 

14 

Zr/Nb  Alloy 

5 

UF6  Nozzle  and  MHD  Channel 

15 

Be 
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16 

A1203  or 
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10  UF6  Core 


Fig.  5.2  Two  dimensional  model  for 
SQUID  calculations. 
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SQUID-360  DIFFUSION  CALCULATIONS 


3  CM  AL203 
6  CM  A L 2 0 3 
3  CM  BEO 
6  CM  BEO 


NOZZLE  AND  MHD  CHANNEL  REGIONS  (  CM 


Fig.  5.4  Integrated  power  for  each  region. 
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time  ns 


Fig.  5.5  Gas  density  5.10el7  molecules/cm3 . 

Temperatures  in  eV.  nf=l . 25. lOelO 
fissions/s. 


time  p.8 


Fig.  5.6  Gas  density  3.10el8  molecules/cm3. 
Other  conditions  as  in  Fig.  5.5. 


Fig.  5.7  Gas  density  10e20  molecules/cm3. 
Note  the  ns  scale. 


neutral  density  (10  /cm  ) 

Fig.  5.8  Stationary  values  of  Tp,  Ts=Ti,  R,  and 
vs.  gas  density. 


-3 


Fig.  5.9  Equilibrium  chemical  composition  of 
UF6  at  10  atm.  (from  Ref. 8). 
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Fig.  5.11  Stationary  ionization  rates  vs.  position. 
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6  SUMMARY  AND  CONCLUSIONS 

This  report  presents  the  results  o-f  work  performed  under  DOE 
contract  number  DE-AC22-88PC88651 .  The  overall  objective  of  this 
project  was  to  investigate  the  use  of  a  gas-core  nuclear  reactor 
directly  coupled  to  an  MHD  expansion  channel  in  an  SDIO  mission 
application.  In  the  proposed  concept,  working  gas  could  be  either 
uranium  hexafluoride  or  a  mixture  of  uranium  hexafluoride  and 
helium,  while  reactor  control  is  performed  by  changing  the 
conf iguration  of  a  reflector  located  outside  the  pressure  vessel. 
Work  was  focused  on  the  basic  physical  mechanisms  of  neutronics, 
reactivity  control ,  ionization,  and  electrical  conductivity  of 
the  concept.  Specific  goals  of  the  project  were: 

-  To  evaluate  the  performance  of  the  external  reflector. 

-  To  determine  the  influence  of  additional,  fission-induced 
mechanisms  on  the  ionization  levels  of  the  working  fluid. 

-  To  estimate  the  electrical  conductivity  levels  that  can  be 
attained  in  the  MHD  channel,  taking  into  account  the  effects  of 
F-  or  other  negative  ions. 

Task  1.,  Reflector  Modeling,  provided  a  reactor  configuration 
that  was  used  as  a  reference  for  subsequent,  detailed  simulation. 
This  configuration  is  featured  by  a  cylindrical  reactor  cavity, 
surrounded  by  a  couple  of  reflecting  (and  refractory)  shells 
located  inside  the  pressure  vessel  and  by  a  movable  sliding 
reflector,  located  outside  the  pressure  vessel.  Two  conical 
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segments,  at  each  end  of  the  cylindrical  core,  complete  the 
reactor.  A  simple  axial-inlet  flow  pattern  has  been  adopted  in 
the  analysis. 

A  number  of  important  features  of  the  proposed  gaseous  reactor 
were  investigated.  The  nuclear  analysis  addressed  the  following 
items: 

a)  materials  selection  for  the  fixed  reflectors  inside  the 
reactor  pressure  vessel,  and  for  the  movable  reflector; 

b)  reactivity  control  performance  of  the  movable  reflector; 

c)  power  density  distributions  inside  the  core  accounting  for  the 
interaction  with  local  fuel  density  and  temperature,  as 
resulting  from  thermal-fluid  analysis. 

d)  reactivity  requirements. 

Task  2. ,  Working  Fluid  Characteristics,  addressed  the  study  of 
the  ionisation  levels  in  the  working  gas,  identifying  a  mechanism 
to  convert  the  energy  released  by  the  fission  reactions  into 
charged  particles,  and  providing  a  system  of  rate  equations  for 
the  numerical  estimates  performed  in  Task  3.. 

Task  3.,  Electrical  Conductivity  Effects,  provided  the  evaluation 
of  the  conductivity  levels  attainable  in  the  generator. 

To  enable  this  evaluation,  the  nuclear  analysis  was  extended  to 
include  the  estimate  of  fission  and  gamma  power  densities 
inside  the  nozzle  and  the  MHD  channel,  taking  into  account  the 
effects  of  the  wall  material  and  thickness.  The  rate  equations 
were  solved  for  a  number  of  typical  gas  conditions  and  duct 
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lengths.  The  influence  of  dissociation  and  ion  attachment,  not 

included  in  the  simplified  numerical  approach,  was  also 

discussed. 

Two  major  conclusions  can  be  drawn  from  the  accomplishment  of 

this  project: 

-  A  reactor  configuration  has  been  obtained  which  avoids  to 
decouple  the  neutronics  of  the  cavity  and  its  internal 
reflector  from  the  external  shell,  thus  enabling  the  latter  to 
be  used  -for  reactivity  control. 

In  particular,  the  results  show  that  the  movable  reflector  is 
able  to  make  the  reactor  subcritical  by  1000  pern  at  least,  to 
control  the  reactor  from  CZP  to  HFP  condition,  accounting  for  a 
reactivity  change  up  to  3300  pem  associated  to  depletion. 

-  The  conditions  expected  in  the  MHD  channel  show  that  the 
ionization  and  conductivity  levels  can  be  amplified  through 
■secondary'  processes  induced  by  the  'primary'  hot  electrons 
which  are  directly  generated  in  the  fission  reactions. 

The  numerical  calculations  suggest  that  conductivity  levels  in 
the  range  10-100  <»m>-«  are  consistent  with  the  characteristics 

of  short  nozzle,  high  power  density  generators.  In  fact,  the 
nozzle  length  turns  out  to  be  of  critical  importance  because  of 
the  sharp  decrease  of  residual  fission  and  gamma  densities  as 
a  function  of  the  distance  from  the  reactor  outlet. 

This  conclusion  remains  true  even  if  the  effects  of  dissociation 
and  attachment  are  included  in  the  numerical  model. 
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Furthermore,  a  preliminary  evaluation  of  the  influence  of  wall 
material  and  thickness  shows  the  potential  for  the  fission  and 
gamma  densities  in  the  duct  to  be  increased  of  one  order  of 
magnitude  at  least,  in  comparison  with  the  values  used  for  the 
conductivity  calculations,  if  a  material  with  favourable  nuclear 
properties  is  adopted. 
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